Pseudo-classical theory for fidelity of nearly resonant quantum rotors 
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Using a semiclassical ansatz we analytically predict for the fidelity of 5-kicked rotors the occurrence of 
revivals and the disappearance of intermediate revival peaks arising from the breaking of a symmetry 
in the initial conditions. A numerical verification of the predicted effects is given and experimental 
ramifications are discussed. 
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Besides entanglement in multipartite systems, it is the 
evolution of phases and the superposition principle which 
distinguishes a quantum from a classical system. Phase 
evolutions can be monitored in many ways, e.g., by corre- 
lation functions 1] . A quantity which has gained interest 
in the last decade is fidelity (2j , defined as the overlap of 
two wavefunctions subjected to slightly different tempo- 
ral evolutions. The temporal evolution of this quantum 
fidelity crucially depends on evolving relative phases. For 
many-particle systems, fidelity can be viewed as a Hilbert 
space measure to study quantum phase transitions Q 
and the regular-to-chaotic transition in complex quan- 
tum systems For single-particle evolutions fidelity 
was measured in electromagnetic wave j&j and matter 
wave @ billiards, and with two different methods for pe- 
riodically kicked cold atoms @, [1] ■ 

The latter system is a realization of the quan- 
tum kicked rotor (QKR), the standard model for low- 
dimensional quantum chaos and the occurrence of dy- 
namical localization Q. Great interest in the QKR has 
reemerged in the study of its quantum resonant motion 



ccmcrgcd m the study ot its quantum resonant motion 
^PJpJlljB,[il,[il] and related accelerator modes 
These two regimes are far from the 
classical limit of the QKR and, therefore, governed by 
distinct quantum effects. Nevertheless, close to quan- 
tum resonance the system can be described (pseudo- 
)classically with a new Planck's constant, which is the 
detuning from the exact resonant value of the kicking 
period [13, [H, [H|- For the quantum resonances, the un- 
derlying pseudo-classical model is completely integrable 
and corresponds in good appro ximation to the dynamics 
of a classical pendulum [Hj, |2l[ . 

In this paper we apply well-known semiclassical meth- 
ods to describe the behavior of fidelity close to the lowest- 
order quantum resonances of the QKR. We extend previ- 
ous analytical results at exact resonance to a broader 
parameter regime, recently measured in experiments per- 
formed by Wu and co-workers B. The behavior of clas- 
sical [22| and quantum fidelity [23 . [24| , in the case when 
classical motion is integrable, has mainly been addressed 
numerically so far, while our approach is both numeri- 
cal and analytical. Also, the recurrences of fidelity found 
in [23| for the near-integrable regime of the kicked ro- 
tor are just predicted for perturbative variations around 



small kicking strengths. Our results are more general, 
allowing, e.g., for strong changes of the fidelity param- 
eter as long as the motion remains nearly resonant. As 
expected, in the nearly-resonant regime, the temporal be- 
havior of fidelity follows the behavior at exact resonance 
the longer, the smaller the detuning from resonance. In- 
deed, we show that the exactly resonant result, predicted 
in by quantum calculations, is retrieved by pseudo- 
classical analysis. At large times, however, the exactly 
resonant fidelity and the nearly resonant one differ, as 
the latter displays recurrent revivals, while the former 
steadily decays. Such revivals are approximately peri- 
odic. Their period depends on the detuning from res- 
onance and diverges as exact resonance is approached, 
so this noteworthy phenomenon is unrelated to quantum 
resonant dynamics. On the other hand, it is quite unex- 
pected on classical grounds because the system is chaotic 
in the proper classical limit. Revivals of fidelity are thus 
a quantum effect and yet are explained by a (pseudo- 
classical analysis that relates them to periodic motion 
inside pseudo-classical resonant islands. Experimental 
possibilities to verify our predictions are discussed at the 
end of the paper. 

The dynamics of kicked atoms moving along a line in 
position space is described, in dimensionless units, by the 
Hamiltonian [13, HI: 



H(t) 



fccos(x) Kt-t'), (!) 



where x is the position coordinate and p its conjugate 
momentum. We use units in which ft = 1 so the param- 
eter r plays the role of an effective Planck's constant; t 
is a continuous time variable, and t' is an integer which 
counts the number of kicks. The evolution of the atomic 
wave function ip(x) from immediately after one kick to 
immediately after the next is ruled by the one-period 
Floquet operator C/fe = exp(— ik cos(z)) exp(— irp 2 /2). 
Fidelity of the quantum evolution of a state ip with re- 
spect to a change of the parameter k from a value k± to 
a value ki is the function of time t which for all integer t 
is defined by: 



F{k u k 2 ,t) = k&imK^I 



(2) 
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Periodicity in space of the kicking potential enforces 
conservation of quasi-momentum (3, which is just the 
fractional part of p thanks to ft = 1. The atomic 
wave function decomposes into Bloch waves [HI, 
which are eigenfunctions of quasi-momentum, ip(x) — 
J^d/3 e lf3x y/plj) #/j(0), where 6 = x mod (2tt) and 
the factor p((3) is introduced in order to normalize ^ p 
(it weights the initial population in the Brillouin zone 
of width one in our units). The dynamics at any fixed 
value of (3 is formally that of a rotor on a circle, param- 
eterized by the angle coordinate 9 and described by the 
wave function ^fp. The Floquet propagator for the rotor 
is given by Up t k = exp(— ikcos(9)) exp(— ir{N + (3) 2 /2), 
where 



-i-jg- Fidelity @ may then be written 



F(k 1) k 2 ,t) = 



dl3 P {p){Ul k J>p\Ul k ^t 



(3) 



so it results from averaging the scalar product under the 
integral sign over (3 with the weight p{(3)- Note that the 
rotor's fidelity is the squared modulus of this quantity, so 
the fidelity ([2]) of atomic evolution does not coincide with 
the /3-average of the rotors' fidelities, c.f. [TTj . Whenever 
r = 2ir£ (£ integer), the evolution is explicitly solvable 
[To| and in particular the rotor's fidelity is determined 
by [Hi 



= Jli\Wt\6k) 



(4) 



where Jo is the Bessel function of 1st kind and order 0, 
8k = k 2 ~ ki and \W t \ = \ sm{nt£(l3 ~ \) csc(n£((3 - |)|. 
If 2/3 — 1 is an integer then a so-called QKR resonance 
occurs and eq. ^ decays in time proportional to t . 
When t is close to a resonant value: r = 2tt£ + e, the 
quantum rotor dynamics may be viewed as the formal 
quantization of t he p seudo-classical dynamics, defined by 
the map 0, EE Hf : 

It+i = It + fcsin(0 t+ i) , 

t+ i = 9 t + I t + ttC + t(3 mod (2tt) , (5) 

using e as the Planck's constant, I — eJV, and k — ek. It 
is thus possible to investigate the quantum fidelity in the 
limit of small e by means of standard methods of semi- 
classical approximation. In the limit e — ► 0, the physical 
parameter k is fixed, so the pseudo-classical parameter 
k — > 0. As a consequence, for sufficiently small e, the 
pseudo-classical dynamics ([5]) is in the quasi-integrable 
regime, even in cases when the classical kicked rotor dy- 
namics is fully chaotic. It is dominated by the resonant 
islands at I res = (2m + £)tt — r/3, with m integer. As 
we consider initial atomic states with a narrow distri- 
bution of momenta near p — 0, we may restrict our- 
selves to a portion of the pseudo-classical phase space 
that includes the one island which is located astride 
1 = 0. We assume I = 1 for simplicity. The pseudo- 
classical dynamics inside the resonant island is ruled, 



in continuous time, by the pendulum Hamiltonian [26[ 
H(6,I,k) = \{I + f3) 2 + fccos(0), where /? = T(/3-i). 
We choose 1^(0) = (2^Y 1 I 2 



1 



2tt 



E 



06 
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so 

-1/2 



,-§<M0,t)-i$^ 



(6) 



where e > is assumed with no limitation of generality. 
The sum is over all trajectories (labeled by the index s) 
which start with I = at time t = and reach position 
9 at time t. 9' = 9' s are their initial positions, and the 
function whose derivative is taken in the pre-factor yields 
9 at time t as a function of position 9' at time 0, given 
that the initial momentum I' = 0. Finally, the function 
$ s {9, t) = S(6,6' s ,t) is the action of the s-th trajec- 
tory and v s is the Morse-Maslov index 27[ ■ We restrict 
ourselves to librational motion inside the stable island. 
The frequency of this motion decreases from u = Vk at 
the island center to lu = at the separatrix. For times 

less than the minimal half-period n/Vk, there is a single 
trajectory in eq. ([5]). Furthermore, 



d'i> B (9,t) _ ( as(e,e',t) ee'(e,t) 
at ~ \ ae 1 dt 



+ 



dS(s.e',t) 
at 



e>=8' 



( dS(6,8',t) 
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= -H( 



,o) =-£ 



fccos(^).(7) 



For t fixed and e — 
equation, so $(0, i) 



we use 9'(9,t) 
r f t 
Jo 



9 - fit in this 
k /J dt' cos(0 - 0t) = 



-\l3 2 t + 2|sin(fi)cos(0 - ^). Replacing all this in 
eq. ([3]), we find for the rotor's fidelity in the limit when 
e — > at constant t: 
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2tt J 
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d Q e iB{0,t)cos{e-^) 



= Jo(B(M) . (8) 
Since B(j3,t) w \W f \ in 



2f sin(f) 



where B(f3, t) — 2 

eq. ((H) for r = 2ir£ and (3 w i, we see that the pseudo- 
classical approximation along with the pendulum approx- 
imation well reproduce the exact quantum calculation (T?|) 
when e — > at fixed t. In the final step of integrating over 
quasi-momenta to find the fidelity for atoms (as distinct 
from fidelity for rotors), the pseudo-classical approxima- 
tion plays no role since the particle's dynamics, unlike the 
rotor's, does not turn pseudo-classical in the limit e — > 
[l8| ■ Replacing (TJJ in eq. ([3]) and computing the integral 
with a uniform distribution of (3 in [0, 1) shows that the 
complete fidelity ([3]) saturates to a non-zero value in the 
course of time [llj |. 

Next we address the asymptotic regime where e — * 
and te 1 / 2 ~ const. To this end, the exact solution of 
the pendulum dynamics is needed in order to compute 
actions; however, some major features of fidelity are ac- 
cessible by exploiting the harmonic approximation of the 
pendulum Hamiltonian. We replace the pendulum by the 
quadratic Hamiltonian H(I, 9) = i(/+/3) 2 + i j-0 2 , where 



3 



and a shift of 9 by tt is understood. Except at ex- 
act multiples of the period, there is one harmonic oscilla- 
tor trajectory in the sum in ([5]); moreover, Maslov indices 
do not depend on the trajectory. Straightforward calcu- 
lations yield 9'(8,t) — sec(a>i)(6l — fluj^ 1 sin(W)) and 
t) = (39 (sec(wt) - 1) - (cj- 1 /? 2 + Lu9 2 )tan{ujt)/2, 
and so: 



j\(t) 



2lT\J\ COs(wit) C0S(CJ2^)I " — 7T 



d6» e ^ 



^-{A(t)9 2 +C(t)/3 2 -2,3eS(t)} 



(9) 



where A(t) = u> 2 ta,n(uj 2 t)—uJi tan(wit), B(t) = sec(w2i) — 
sec(wii) and C(t) — lj^ 1 t&n(uj 2 t) — wf 1 tan(wii). A(i) 
is a phase factor accumulated by the Maslov indices and 
it just depends on time, rendering it irrelevant for our 
present purposes. We next insert eq. © in eq. (J3]) and 
choose for p(f3) a uniform distribution in some interval 
[| — b, | + 6Jj with < b < 1/2. It is necessary to as- 
sume that b is smaller than the halfwidth of the pseudo- 
classical resonant island, because the harmonic approx- 
imation we have used is valid only inside that island. 
Then 



F(k u k 2 ,t) 



1 



l6ir 2 b 2 T 2 1 cos (tuit) cos (u> 2 t) \ 



d8 e~~c 



A 1 (9,e,t) 



d/3e~ 



-A 2 (0,S,e,t) 



-rb 



(10) 



and 



where Ai(0,e,t) = (A(t) - B 2 (t)C{t)- r ] 
A 2 (P,9,e,t) = ((3^/C(t) - B{t)C{t)- 1 / 2 9) 2 . As A 2 ~ 
e -!/ 2 j n the limit when e — ► and t\J~i ~const., the lim- 
its in the ^-integral in (JTUJ) may be taken to ±oo: 



d/3 e~2 



£A 2 (/3,0,e,t) 



(27T) 



1/2,1/2 



C{t) 



-l/2 e -i7r/4 



Due to this approximation, ()11|) below is valid in the 
asymptotic regime where e is small compared to b 2 . The 
remaining ^-integral is dealt with similarly, because the 
pre-factor of 9 2 in Ai is ~ e~i/ 2 . Thus finally 

e 2 

F(k u k 2 ,t) ~ 167T 2 b 2\ c ^ A ^ _ s(t)2|| cog^t) C0S (u;2i)| 



87r 2 6 2 |4o;iaj2 — u\ cos(ui-t) — lo 2 _ cos(w+i)| 



(11) 



where lu± — lo\ ± lo 2 . Singularities of this expression are 
artifacts of the approximations used in evaluating the in- 
tegrals in (|10p , which indeed break down when the divisor 
in (|11[) is small compared to e. However, they account for 
the periodic "revivals" that are observed in the fidelity 
at large times, with the beating period T\ 2 = 2tt/\uj^\ 
(Fig. 2 (a)). With a quite narrow distribution of (3, how- 
ever, fidelity is at long times dominated by the "resonant" 
rotors {(3 = or j3 = 1/2 respectively), and then revivals 
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FIG. 1: (color online) Fidelity as predicted by eq. (|12p - be- 
cause of the singularities of the analytical formula the curve 
is folded with normalized Gaussians with a standard devia- 
tion of t ~ 6 kicks (solid black line) - and numerical data 
(grey/green curve), for fci = 0.87T, hi = Q.6n and detuning 
e = 0.01 from r — e = 2n. In the inset, the non-smoothed 
result (1121) is shown. 



occur with the period T 12 /2 (Fig. 1). 
purely resonant (3, eq. (fTU|) yields: 



Indeed, with the 



F(k 1 ,k 2 ,t) = F res (k 1 ,k 2 ,t) 
1 



2tt \uj 2 cos(wii) sm(uj 2 t) — cos(u> 2 t) sin(wii) | 



, (12) 



which has singularities in time with the mentioned peri- 
odicity of T12/2. This behavior of resonant rotors has a 
simple qualitative explanation. As the initial state of the 
rotor corresponds to momentum I — 0, at that value of 
quasi-momentum (fi — \) the stationary-phase trajecto- 
ries of the two harmonic oscillators, which were started 
at I = Q, exactly return to / = whenever time is a 
multiple of the half-period T± 2 /2, and so fully contribute 
to fidelity, in spite of their angles being different by tt in 
the case of odd multiples. At j3 ^ this symmetry is 
lost. Comparing numerical data (obtained by repeated 
application of the Floquet operator to the initial wave- 
function) with the analytical predictions we find excellent 
agreement. We observe the expected peak structure of 
the revivals in Fig. [T] and the loss of intermediate revival 
peaks at T\ 2 j2 in Fig.^a). The time scale on which the 
revivals occur is proportional to e" 1 / 2 and of crucial im- 
pact to experimental measurements: conservation of co- 
herence has been shown for up to 150 kicks (see [19() with 
cold atoms, making an observation of the revivals for rea- 
sonable e < 0.01 possible. Earlier realizations of the QKR 
were implemented using cold atoms @, HI, 13 with broad 
distributions in quasi-momentum. Nowadays, much bet- 
ter control of quasi-momentum is provided by using Bose- 
Einstein condensates (see [13j, LL6|), which allows for a 
restriction in (3 up to 0.2 % (as achieved in [lj|) of the 
Brillouin zone. This would allow to verify our results by 
conveniently reducing the intervals in quasi-momentum 
and thus retracing the revivals with period T 12 /2 to the 
exactly resonant and the revivals with period T\ 2 to the 
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FIG. 2: (color online) Same as in Fig. [T] for an ensemble of 
5000 equidistantly chosen rotors (solid grey/green lines) with 
a width of (a) A/3 = 0.05 (or A/3 » 0.31) around the reso- 
nant value, covering half the width of the resonance island in 
the phase space induced by ©, and (b) A/3 = 1, covering 
the full phase space, compared with the smoothed (see cap- 
tion of Fig. [TJ version of eq. (solid black lines). In (a) 
the intermediate revival peaks observed in Fig. Q] disappear 
as predicted by (|11[>. The dashed line in (a) reproduces the 
smoothed analytic formula from Fig. [T] For /3 distributed over 
the full Brillouin zone in (b), the revivals are barely visible 
since the average includes many nonresonant rotors perform- 
ing rotational motion in phase space, which is not decribed 
by our theory valid just for the librational island motion. 



near-resonant rotors. There exists an interesting second 
possibility to measure the transition from eq. (|12[) to 
eq. (fTT|) with just cold atoms, since the (3 we use scales 
with the kicking period, i.e. f3 — r(/3 — 1/2). Due to this 
scaling, the limit r — ► (automatically implying also 
e — > 0, c.f. [14J) permits a measurement of eq. (|12|). even 
with an ensemble of cold atoms whose quasi-momenta 
occupy the full Brillouin zone. Also the momentum se- 
lective interferometric measurements of fidelity Q allow 
to select narrow intervals of quasi-momenta, and hence 
would permit to check our predictions experimentally. 

To summarize, we predict fidelity revivals in the QKR 
close to quantum resonance using a semiclassical ansatz. 
Our results are supported by numerical data showing 
the same characteristic revival peaks. Every second 
peak vanishes once the symmetry of the initial quasi- 
momentum distribution on the resonance island is bro- 
ken. This makes for a surprising transition that could be 
measured with both cold and ultracold atoms owing to 
the scaling of (3 or the use of momentum selective meth- 
ods, as described in the previous paragraph. 
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